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The processes of six-quark production with one bb pair are studied by means of a 
^ ' complete tree-level electroweak calculation. The top-quark signal is examined: the 

• importance of electroweak backgrounds, of the order of 10% above the tt threshold 

and of about 30% of the purely electroweak signal at threshold, is further stressed 
by studying the dependence of the cross-section at threshold on the Higgs mass 
in the range between 100 GeV and 185 GeV, and finding variations of the order 
of 10%. In the study of some event-shape variables, a strong effect of initial-state 
radiation is found, in particular for the thrust distribution, which is studied for 
several centre-of-mass energies at the TeV scale. The effectiveness of cuts on the 
thrust for isolating QCD backgrounds, as pointed out by some authors, is confirmed 
also in the presence of electroweak backgrounds and initial-state radiation. 
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1 Introduction 



Many important signals to be studied at future high energy coUiders, NLC [1] 
and LHC [2], will have a large number of particles in the final state. In particu- 
lar, the processes with six fcrmions in the final state will be of great importance 
for several tests of the Standard Model, such as the studies of top-quark and 
electroweak gauge bosons, as well as the search for an intermediate-mass Higgs 
boson. Such processes are already of great interest at the Tevatron collider [3] 
in connection with top-quark physics. Theoretical studies of six-fermion (6/) 
processes by means of complete tree-level calculations have only very recently 
appeared in the literature [4-9], where top-quark physics, Higgs physics and 
WWZ production have been addressed. 

All these studies clearly demonstrate that complete calculations are important 
for a precise determination of the cross-sections, and for the development of 
reliable event generators, whenever accurate evaluations of interference, off- 
shellness and background effects, as well as spin correlations, are important. 

In view of the precision measurements of the top-quark properties we have 
to analyse, among the 6/ signatures, the ones containing a bb pair and two 
charged currents, as the top-quark decays almost exclusively into a W bo- 
son and a b quark. Semi leptonic signatures have already been considered in 
refs. [4,5,10]. It is then of great interest to carefully evaluate the size of the 
totally hadronic, six-quark (6g) contributions to integrated cross-sections and 
distributions as well as to determine their phcnomenological features. The aim 
of the present study is to make a first quantitative analysis in this context, 
for what concerns the full set of electroweak contributions to a class of e+e~ 
annihilation processes related to top-quark physics. Special emphasis will be 
given to the determination and to the analysis of the topology of the events 
considered, so as to characterize them, as far as possible, against the QCD 
backgrounds. 

Looking at an experimental situation where the 6-tagging technique can be 
applied, it is meaningful to distinguish the 6q final states containing one bb 
pair, of the form bbqq'q"q"', from those containing two or three bb pairs, re- 
spectively of the form bbbbqq and bbbbbb. The last two kinds of processes are 
not relevant to top-quark production, as they contain no charged currents. 

In the present study the signatures with one bb pair are considered and the full 
set of purely electroweak contributions is taken into account. These processes 

can be further divided into three subsets (although in realistic predictions 
they cannot be treated separately), which are shown in Table 1: concerning 
the quark fiavours other than b, only charged currents are involved in the first 
subset, both charged and neutral currents in the second one, and only neutral 
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bbddss 



Table 1 

Six-quark final states with one bb pair. The notations CC (charged currents) and 
NC (neutral currents) refer to the currents formed by the quark flavours other than 




Fig. 1. The two Feynman diagrams with tt production. 

currents in the third one. The total number of tree-level Feynman diagrams 
involved in the complete electroweak calculation amounts to several hundreds. 
Such a complexity is unavoidable, as will be shown, if an accuracy of 1% is to 
be reached. The diagrams with top-quark production, which will be referred 
to as signal diagrams, are shown in Fig. 1. They contribute to the processes 
in the first two columns of Table 1, but not to those in the third. 

All the processes receive contributions from diagrams of Higgs production, 
of which the leading ones are illustrated in Fig. 2. The relevance of such 
contributions depends on the Higgs mass and on the centre-of-mass (cm.) 
energy: the dominant decay mode is H ^ bb for low Higgs masses {rriH < 130- 
140 GeV) and H VV (V = W, Z) for high Higgs masses. The predictions 
can thus be expected to depend on the Higgs mass. 

As is well known, the behaviour of the cross-section near the threshold for ti 
production is characterized by strong interaction effects that give a sizeable 
modification with respect to the purely electroweak prediction. Such effects 
are treated in the literature [11,12], and are not included in the calculations 
presented in this paper. Q Results at energies around the threshold are shown, 

^ Theoretical calculations of radiative corrections to ti production are also present 
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Fig. 2. Leading Higgs contributions, divided into charged-current terms, on the left, 
and neutral-current terms, on the right (where charged and neutral currents are 
referred to quark flavours other than b). Other diagrams, obtained by attaching the 
Higgs line to the other possible gauge-boson line in the diagrams (d), (e), (/), and 
(g), are understood. 

so as to give a thorough analysis of the electroweak contribution. Some of the 
QCD backgrounds to the signatures considered in the present study have been 
evaluated in ref. [13], and their topology has been studied by means of event- 
shape variables. One of the objectives of the present work is to characterize 
the topology of the complete electroweak contributions in order to help finding 
appropriate selection criteria to reduce as far as possible the QCD backgrounds 
studied in ref. [13]. The analysis performed in the present work, together with 
the other studies in the literature so far, should give a complete picture of 
electroweak contributions to 6/ processes relevant to top-quark physics at 
NLC. 



in the literature, as recently reviewed in ref. [12]. 
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The paper is organized as follows: in Section 2 the computing procedure is 
briefly described; in Section 3 the numerical results, including integrated cross- 
sections and various distributions, are presented and discussed; Section 4 is 
devoted to our conclusions. 



2 Calculation 



The numerical results have been obtained by means of a procedure analogous 
to the one adopted in ref. [8], where the interested reader can find some tech- 
nical details that will be omitted here. The computer program already used 
in ref. [8], which is based on ALPHA [14] for the matrix element calculation 
and on an evolution of HIGGSPV/WWGENPV [15,16] for the Monte Carlo 
integration and event generation, has been adapted, in the multichannel im- 
portance sampling, to include some new diagram topologies, such as those in 
Figs. 1 and 2. 

The cross-section is calculated according to the formula 

cr = [ dzidz2DBs{zi,Z2;s) [ dxidx2D{xi, s)D{x2, s)d[PS] f^r,i > (1) 



where initial-state radiation {ISR) [17] and beamstrahlung (BS) [18] are in- 
cluded by means of the structure functions D{x, s) and Dbs{x-i s), respectively; 
da/d[PS] is the differential cross-section at the partonic level, and d[PS\ is 
the six-body phase-space measure. The program may be used to generate un- 
weighted events as well. 

The input parameters are G^, M^, M^, the top-quark mass rut — 175 GeV, 
and the 6-quark mass = 4.3 GeV; all the other fermions are treated as 
massless. The widths of the W and bosons and of the top-quark and all 
the couplings are calculated at tree level. The Higgs-boson width includes the 
h — > jiji, TT, cc, hh, the h ^ gg [19] and the two- vector-boson [20] channels. The 
CKM matrix used is exactly diagonal. The propagators of unstable particles 
have denominators of the form — + irM with fixed widths. The validity 
of this choice for minimizing possible gauge violations has been discussed in 
ref. [8], where the final states qql^l'vi' were considered. In that paper, for 
the SU{2) invariance, the fudge-factor method has been used to check the 
numerical results: apart from the well-known problems of the fudge scheme, 
i. e. the mistreatment of non-resonant diagrams close to the resonances, no de- 
viation has been found in the total cross-section up to the numerical accuracy 
considered. In order to check U{1) invariance, the matrix element has been 
calculated with different forms of the photon propagator obtained by varying 
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the gauge parameter; the results were found to be stable up to numerical pre- 
cision. The same analysis carries on to the present study and gauge-violation 
effects are estimated to be numerically negligible. 

The colour algebra, not implemented in the version of ALPHA that has been 
employed here, has been performed by summing the different processes with 
proper weights. As an example of this, the process e^e~ hhudud may be con- 
sidered: the colour amplitude, in the case of purely electroweak contributions, 
can be written in the form 

A = {aiSi^ir^Si^i^ -\- 0251113^2214) ^jk -I (2) 

where the colour indices 11,12, Hi U, j and k refer to the u, d, u, d, b and b quarks 
respectively. The squared modulus summed over colours is then 

j2 i^r = i«ir + («i«2 + «i«2) + i«2r . (s) 

col 

The amplitude given by ALPHA is instead 

A^ai + a2 . (4) 

Thus one cannot use an overall factor to obtain eq. (3) from eq. (4). In order 
to disentangle the various terms in eq. (3), it is useful to notice that, with the 
quark masses adopted here and with a diagonal CKM matrix, the first term 
in the right-hand side of eq. (2) is equal to the amplitude A' of the process 
e+e~ — > bbudcs, and the second term is equal to the amplitude A" of the 
process e+e~ — > bbusus. Similarly, the two colourless amplitudes A' and A" 
of these processes arc equal to the first and to the second term, respectively, 
in the right-hand side of eq. (4). Thus the following relation is valid: 

J2 (l^r + l^f + = (l^r + - l){\Af + \A'f)) ■ (5) 

col 

Other situations are treated in a similar way, and the correct colour weights 
are thus obtained in the sum over the whole class of processes considered. 



3 Numerical results and discussion 

In this section the numerical results, including both integrated cross-sections 
and distributions, are shown. In all the calculations the invariant masses of 
the bb pair and of all the pairs of quarks other than b and b are required 
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to be greater than 10 GeV. The results presented below arc obtained, unless 
otherwise stated, by summing over all the processes listed in Table 1. 

3.1 Integrated cross-sections 

As a first step, the total cross-section, resulting from all the tree-level diagrams 
for the processes in Table 1, has been calculated in the Born approximation at 
energies from 340 to 800 GeV. Two values of Higgs mass have been considered, 
rrifj = 100, 185 GcV, so as to study the dependence of the results on in 
the intermediate range. The numerical errors are always below 1% and in 
particular above the tt threshold they are kept at 0.2 — 0.3% level. 

In Fig. 3 the full cross-section for mn = 185 GeV is compared with the signal, 
defined as the contribution of the two diagrams of ti production of Fig. 1, 
summed over the four processes to which they contribute (see the first two 
columns of Table 1); the signal is shown both in the Born approximation and 
with ISR switched on. 

The difference between the full and the signal curve is dominated by Higgs- 
strahlung contributions (diagrams (a) and (6) of Fig. 2) at low energy, while 
other backgrounds are important at high energy, coming from all the processes 
in Table 1 and, for a little amount, from the interference of the signal diagrams 
with the other contributions in the chargcd-current and mixed processes. The 
electroweak background effects, which are in the range 5 — 10% above the 
threshold, amount to around 30% at threshold; they are much greater below 
the threshold, where the signal is suppressed with respect to the background 
(the ratio background/signal is of 2.5 at 340 GeV). 

The radiative effects strongly suppress the cross-section in the low-energy 
region, where it grows rapidly, as they reduce the effective cm. energy; with 
increasing energy, the curve with ISR comes to cross the one in the Born 
approximation, as a consequence of the onset of the opposite behaviour of the 
Born term, which, above the threshold, starts decreasing. It can be observed 
that at 500 GeV the enhancement due to the background is of the same order 
as the lowering given by the ISR. 

In Fig. 4 the signal cross-section without kinematical cuts is plotted together 
with the cross-section in the narrow- width approximation (NWA). The latter 
is calculated as the product of the cross-section for e+e~ — > tt, and of the 
branching ratios of the decays W qq', assuming the branching ratio of 
t — * Wb to be exactly unity. The difference between the two calculations is 
about 15% in the region near the threshold, and it decreases, as expected, 
with increasing cm. energy: at 500 GeV it is 3%, while at 800 GeV it is less 
than 1%. These results give a measure of the off-shellness effects connected 
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Fig. 3. Full six-quark electroweak cross-section (solid line) and tt signal (dashed line) 
in the Born approximation, and tt signal with initial-state radiation (dash-dotted 
line), as a function of the cm. energy. 




700 750 BOO 

V s (GeV) 



Fig. 4. Signal cross-section (solid line) without cuts compared with the narrow- width 
approximation (dashed line), as a function of the cm. energy. 

with the top-quark and VF-boson widths. 

The cross-sections for the two values of the Higgs mass, rriH = 100 and 185 
GeV, have been found to differ by less than 1% at energies above the threshold 
region, while at lower energies, differences of up to 20 — 30% occur. This is 
due to the fact that the signal at low energy is not large enough to hide the 
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Vs = 350 GeV 




100 110 120 130 140 150 160 170 180 

itIh (GeV) 

Fig. 5. Total cross-section as a function of the Higgs mass at the threshold for ti 
production. 

Higgs-mass effects. Moreover, such effects decrease with increasing energy. In 
order to make a detailed study of the dependence on the Higgs mass at low 
energy, the cross-section at the threshold for ti production, ^/s = 350 GeV, 
has been calculated for various Higgs masses in the range from 100 GeV to 
185 GeV. The results are shown in Fig. 5, where the cross-section is plotted 
as a function of the Higgs mass. Variations of the order of 10% can be seen in 
this plot, which shows the importance of complete calculations to keep under 
control the background effects and uncertainties that come from not knowing 
the Higgs mass. 



3.2 Distributions 

Two samples of events have been generated at a cm. energy of 500 GeV and 
with a Higgs mass of 185 GeV. One sample is in the Born approximation, 
while the other includes ISR and BS. The numbers of events, of the order of 
10^, have been determined by assuming a luminosity of 500 fb~^, which is the 
integrated value expected in one year of run. 

In the definition of observable distributions for the class of processes consid- 
ered here, we must take into account the fact that quark flavours other than 
b cannot be identified. As a consequence, two kinds of distributions, labelled 
"exact" and "reconstructed" , are considered in the following: the "exact" dis- 
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tributions are calculated by identifying all the quarks; the "reconstructed" 
distributions are calculated by means of the following algorithm. The mo- 
menta gi, . . . ,g4 of the four quarks other than b and b are first considered 
and, for every pair (gj, qj), the invariant mass TJiij = {qi + qjY is calculated; 
then the two W particles, Wi and W2, are reconstructed as the pairs {qi, qj) 
and (qkjQi) such that the quantity {rriij — Mw\ + Irriki — My[r\ is minimized; 
the top-quark is then determined by taking the combination (6, Wj), {b,Wj), 
which minimizes the quantity Im^Wj — "^t| + I'nT'bWj ~'>^t\, where rht — 175 GeV 
is the nominal top mass. 

The invariant mass of the top-quark is studied in Fig. 6. In the plot (a) a com- 
parison is made between the exact (dashed line) and the reconstructed (solid 
line) distribution in the Born approximation and a good agreement can be 
observed. In order to further check the reconstruction procedure, in particular 
the dependence on the adopted value of rht, some tests have been made by 
taking values in the range 170 GeV < rht < 180 GeV and fitting the resulting 
histograms with Breit-Wigner distributions. The values of the physical top 
mass obtained in the various cases are identical, within the statistical errors. 

The radiative effects are shown in the plot (b) of Fig. 6, for the reconstructed 
distribution. They do not apparently give a substantial modification. In the 
plot (c) the role of background diagrams is studied, by comparing the result of 
the full calculation with the signal alone. The background does not introduce 
any observable distortion. More quantitative results have been obtained by 
means of fits to the histograms with Breit-Wigner distributions. All the his- 
tograms in Fig. 6 give the same value of nit, so that it can be safely concluded 
that electroweak backgrounds, as well as ISR and BS do not give any bias in 
the determination of the physical mass via the direct reconstruction method 
on the scale of precision of 100 MeV. 

The angular distribution of the top-quark with respect to the beam axis is 
a good indicator of the spin nature and of the couplings of the top-quark. 
This variable is illustrated in Fig. 7. As in the case of the invariant mass, the 
exact and reconstructed distributions have been checked to be in very good 
agreement. It should be observed that radiative and background effects are of 
the same magnitude (in particular the former are dominated by the ISR) . The 
shapes of the histograms are in good qualitative agreement with the angular 
distribution predicted by the lowest-order analytic calculation for the process 
of real ti production. 

The most effective way to obtain a separation between the tt signal and the 
QCD backgrounds, as already pointed out by some authors [13,21], is to anal- 
yse event-shape variables, such as thrust [22], sphericity [23], spherocity [24], 
C and D parameters [25], etc. A comparison between pure QCD ((^(ag^ag)) 
six-jet events and the tt signal has been performed in ref. [13] for the thrust 
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Vs = 500 GeV 




172 174 176 178 
rUt (GeV) 

(c) 



Fig. 6. Invariant mass of the top-quark, (a): exact (solid line) and reconstructed 
(dashed line) distribution; (6): distribution in the Born approximation (dashed line) 
and with initial-state radiation and beamstrahlung (solid line); (c): full calculation 
(solid line) and signal (dashed line). 

and sphericity distributions, and other shape variables have been studied in 
the same article for the QCD contributions only. In the present work several 
such variables have been analysed for the electroweak contributions, and the 
effects of the electroweak backgrounds and of ISR and BS have been studied. 
The thrust and C parameter distributions for the process under considera- 
tion are shown in Fig. 8. In the upper row the radiative effects are displayed, 
while in the plots of the lower row the signal is compared with the full result. 
In the radiative case, the distributions are calculated after going to the cm. 
frame. Remarkable effects due to ISR and BS can be seen in these plots and 
in particular in the thrust distribution, where the peak is strongly reduced 
with respect to the Born approximation and the events are shifted towards 
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Vs = 500 GeV 




Fig. 7. Reconstructed angular distribution of the top-quark with respect to the e"*" 
beam axis, (a): results in the Born approximation (dashed histogram) and with ini- 
tial-state radiation and beamstrahlung (solid histogram). (6): full calculation (solid 
line) and signal contribution (dashed line). 

the lower values of T, which correspond to spherical events. It is interesting 
to observe that this phenomenon is of help for the selection of the signal with 
respect to QCD backgrounds. From the plots in the second row, it can be seen 
that the presence of the electroweak backgrounds, although visible, is almost 
negligible for both observables. 

The remarkable change of the thrust distribution after inclusion of radiation 
can be better understood by observing the dependence of this distribution on 
the cm. energy, which is analysed in Fig 9, where four samples of 10000 events 
each, at the energies of 360, 500, 800 and 1500 GeV, are studied. The peak 
structure that is present at 500 GeV is completely lost at 360 GeV, and this 
explains the lowering of the peak at 500 GeV in the presence of ISR and BS, 
as this reduces the available cm. energy. At 800 and 1500 GeV the peak is 
shifted towards the collinear region T" ~ 1, as a consequence of the Lorentz 
boost of the t and i quarks. 

As a conclusion, we can say that, at 500 GeV, in view of the results of the pure 
QCD processes, studied in ref. [13], the thrust variable is the most effective 
in discriminating pure QCD backgrounds, also in the presence of electroweak 
backgrounds and of ISR and BS. At higher energies this separation appears 
to be more and more problematic. 

On the other hand, the backgrounds of 0(«g„a^), given by 2 — * 4 processes 
with subsequent gluon emission from a quark line, should be considered (a 
study of contributions of this class for semi leptonic signatures is made in 
ref. [10], but without an analysis of event-shape variables). A rough estimate 
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Fig. 8. Event-shape variables, (a): thrust distribution in tlic Born approximation 
(dashed histogram) and with initial-state radiation and beamstrahlung (solid his- 
togram); (b): C parameter distribution, as in (o); (c): thrust distribution in the 
Born approximation from the full calculation (solid histogram) and from the signal 
contributions alone (dashed histogram); (c?): C parameter distribution, as in (c). 



of the leading contributions of this kind could be obtained by considering 
a four-fermion process of the form e+e~ W~^W~ — > 4 jets, similar to 
what is done in ref. [21]. A test made by means of the four-fermion program 
WWGENPV has confirmed the results of ref. [21] for the thrust and has led 
to similar conclusions for the C parameter: such processes appear to be well 
separated from the top-quark signal and thus appear to be less dangerous than 
the pure QCD backgrounds. 
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Fig. 9. Thrust distribution in the Born approximation at 360, 500, 800 and 1500 
GeV. 



4 Conclusions 



The production of tt pairs has been studied in processes with six quarks in 
the final state, at the energies of the NLC. The signatures considered contain 

one bb pair, and receive contribTitions from both charged and neutral cur- 
rents. The top-quark signal is present only in the charged-current terms. The 
purely electroweak contributions have been considered and complete tree-level 
calculations have been performed. 

The cross-section has been calculated by means of a computer program already 
used for other phenomenological studies on 6/ processes and adapted here 
to sample the new diagram topologies. The importance of the electroweak 
backgrounds and of the off-shellness effects has been examined. Above the 
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threshold for tt production, the former are of the order of several per cent 
and the latter are at the per cent level. Near the threshold, both effects are 
sizeable and, in particular, a study of the dependence of the cross-section on 
the Higgs mass at threshold shows that variations of the order of 10% occur 
for Higgs masses between 100 GeV and 185 GeV. A complete calculation is 
needed to keep such effects under control and to have a 1% accuracy. 

Some distributions have been studied in a realistic approach, by using a recon- 
struction algorithm for the top-quark that takes into account the impossibility 
of identifying quark flavours other than b. The invariant mass of the top-quark 
has thus been studied and the presence of electroweak background contribu- 
tions as well as the initial-state radiative effects have been found not to affect 
the determination of the mass on the scale of experimental precision expected 
at NLC. 

The angular distribution of the top-quark with respect to the beam axis, which 
is directly related, in the case of real production, to the quantum numbers of 
the top-quark, has been calculated and shown to be in qualitative agreement 
with the expectation suggested by the real production case. 

Finally, some event-shape variables have been studied. At a cm. energy of 
500 GeV the thrust distribution turns out to be the most interesting for the 
aim of discriminating the leading QCD backgrounds, as suggested by other 
authors who discussed the top-quark signal alone. The effects of electroweak 
backgrounds and of ISR and BS have been shown here not to alter these 
conclusions. At higher energies, the Lorcntz boost gives to the event a more 
coUinear shape, so that the separation of QCD backgrounds could become 
more difficult. 

The study presented in this paper has been performed by means of a comput- 
ing program that can equally well deal with semi leptonic signatures and can 
be switched in a straightforward manner to treat polarized scattering. More- 
over, by employing the new version of ALPHA [26] , which embodies also the 
QCD Lagrangian, complete strong and electroweak results could be obtained. 
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